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ABSTRACT 

Previous research showed that most BL Lac objects and some quasars have aligned 
VLBI-core and optical polarizations, although some of the AGNs also showed no ob- 
vious relationship between their VLBI-core and optical polarization angles. This may 
indicate that some AGNs have co-spatial regions of optical and radio emission, while 
others do not. However, another possibility is that some of the VLBI cores had Faraday 
rotations of the order of several tens of thousand of rad/m 2 , which were not properly 
fit using the three-frequency data due to mr ambiguities in the observed polarization 
angles, leading to incorrect subtraction of the effects of the core Faraday rotation, and 
so incorrect intrinsic radio polarization angles xo- With this is mind, we obtained ad- 
ditional 12+15+22+24-1-43 GHz plus optical observations for 8 of 40 AGNs previously 
considered, enabling improved sampling of Faraday effects. Our results indicate that, 
although some VLBI radio cores have comparatively high rotation measures, this alone 
cannot explain the misalignments found between the radio core and optical VLBI po- 
larization angles, Ax = \Xopt — Xo\ - Comparison between A% and (i) the orientation of 
Xo relative to the jet direction, (ii) the degree of polarization of the core, (hi) a depo- 
larization factor, (iv) the core rotation measures and (v) the core magnetic fields 1 pc 
from the jet base do not yield evidence for any correlations between these properties. 
There is, however, some evidence that the maximum observed Ax tends to decrease 
as the core-region magnetic field increases, suggesting that large misalignments in Ax 
could be associated in part with relatively low core magnetic fields. Thus, although the 
overall distribution of Ax for all 40 sources in our sample does show a significant peak 
at Ax ~ 0, it remains unclear what distinguishes these AGN cores from those showing 
appreciable misalignment between optical and VLBI-core polarization position angles. 

Key words: AGN - Faraday Rotation. 



1 INTRODUCTION 

In general, non-thermal synchrotron emission associated 
with the relativistic jets of radio-loud AGNs dominates 
the UV-radio continua of radio-loud active galactic nuclei 
(AGNs). Even so, little correlation between the observed 
emission in widely spaced wavebands has usually been ex- 
pected, even if genuinely simultaneous observations are com- 
pared. This is in part because early attempts to search for 
correlations using integrated single-dish measurements were 
largely unsuccessful (e.g., Rudnick 1978). In addition, opti- 
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cal variations typically occur more rapidly that radio varia- 
tions, suggesting that the optical emission is generated in a 
much smaller volume than the radio emission. In this con- 
text, it was considered natural for there to be little or no 
correlation between the optical and radio emission in AGNs. 

It has only been relatively recently that correlations be- 
tween properties at different wavelengths have been found, 
with several investigations finding evidence for the polar- 
ization properties of sources at different wavebands to be 
correlated. For example, Lister & Smith (2000) carried out 
a joint analysis of VLBA polarization data at 22 and 43 GHz 
and optical polarization data. Despite the fact that the opti- 
cal and VLBA measurements were obtained about a month 
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apart, they found evidence that the optical polarization fol- 
lows the radio structure and tends to lie along the VLBI jet 
direction, with the degree of optical polarization correlated 
with the core polarization and core luminosity at 43 GHz. 

More recently, Gabuzda et al. (2006) carried out 
a simultaneous analysis of optical polarization data and 
15+22+43 GHz VLBI polarization data for 11 BL Lac ob- 
jects and 3C279. Their results showed that, after correc- 
tion for Faraday rotation, nearly all the VLBI core polariza- 
tions were aligned with the corresponding optical polariza- 
tions within 20°. Similar results were obtained by Jorstad 
et. al (2007) at even shorter radio wavelengths. Also, a 
co-rotation of the 7mm and optical polarization angle of 
0420-014 through ~ 80° over about 10 days was found by 
D'Arcangelo et al (2007). 

However, the optical and 15+22+43 GHz observations 
analyzed by Algaba et. al (2010) (hereafter, Paper I) gave 
rise to unexpected results. Based on the work of Gabuzda 
et. al (2006), the initial expectation was a correlation be- 
tween the optical and Faraday-corrected VLBA radio-core 
polarization angles, with the difference between these two 
observed angles, Ax, often being close to zero. As seen in 
Paper I, this is not the case: not only do quasars in general 
not show this trend, but neither do all BL Lacs, in contrast 
with the striking results found by Gabuzda et. al (2006). 

In Paper I we tried to identify some possible reasons 
for the lack of correlation in some sources. In particular, we 
searched for relations between Ax and the degrees of polar- 
ization at the observed radio frequencies, and between Ax 
and the difference between the intrinsic radio-core polariza- 
tion angle and jet direction. No signs of any relation were 
found in either case. 

Paper I discusses possible reasons for this such as the 
existence of internal or high external Faraday Rotation, but, 
unfortunately, the data were not able to unambiguously dis- 
criminate between these possibilities. Although three fre- 
quencies are in principle adequate to fit for the Faraday ro- 
tation and find the intrinsic polarization angle, a variety of 
uncertainties and ambiguities can arise due to factors such 
as the Faraday rotation effects mentioned above, uncertainty 
about optical depth transitions, and beam depolarization. 
We investigate these issues here using five- frequency Very 
Long Baseline Array polarization data for a subset of eight 
of the 40 AGNs. 

2 OBSERVATIONS 
2.1 VLBA 

We obtained a new set of polarization observations on the 
NRAO Very Long Baseline Array (VLBA) at 12.039, 15.383, 
21.775, 23.998 and 43.135 GHz, each with two intermediate- 
frequency bands with a bandwidth of 8 MHz and an aggre- 
gate bitrate of 128 MB/s. These polarization observations 
were carried out in a 24-hour session on November 2, 2008 
(hereafter, epoch D). In all cases, the sources were observed 
in a "snapshot" mode, with 8-10 several-minute scans for 
each object at each frequency spread in time. The result- 
ing coverage in the u-v plane was quite uniform. The data 
reduction and imaging for the radio data were done with 
the NRAO Astronomical Image Processing System (AIPS) 



Table 1. EVPA Corrections 



12GHz 


15 GHz 


22 GHz 


24 GHz 


43 GHz 


[1] 


[2] 


[3] 


[4] 


[5] 


66 + 4 


-63 ±4 


45 + 4 


-85 ±4 


-20 + 4 



using standard techniques. The reference antenna used was 
Los Alamos. Simultaneous solutions for the instrumental po- 
larizations and source polarizations for the compact source 
1954+513 were derived using the AIPS task LPCAL. 

We calibrated the electric vector position angles (EV- 
PAs) using the phases of the VLBA polarization D-terms 
(Gomez et. al 1992). Comparison of the D-terms against 
a set of tabulated values (previously calibrated by other 
means) is a reliable method for calibrating the absolute L- 
R phase offset of the reference antenna in VLBA observa- 
tions. We were able to use two sets of data at the same or 
very similar frequencies, for which the EVPA calibration had 
been performed using integrated polarization measurements 
in the standard way (O'Sullivan & Gabuzda 2009; Reichstein 
& Gabuzda 2011), making it possible to apply this method 
to the data. These other sets of data were obtained on July 
2, 2006 (12.9, 15.4, 22.2, 43.1 GHz) and September 27, 2007 
(12.9, 15.4 GHz). 

To ensure the reliability of this method, we performed 
a series of additional checks. First, we compared the in- 
dependent results given by the two comparison sets of D- 
terms for consistency. Second, we compared the new im- 
ages with previous data (Algaba 2010) for the 8 AGNs in 
common. Third, for some frequencies, we were also able to 
compare the resulting polarization angles with images from 
the MOJAVeQ and Boston Universitjjj databases, as well as 
integrated measurements from the University of Michigaijf] 
database. Fourth, we obtained the EVPAs of several sources 
in the optically thin jet regions and checked for consistency 
with modest Faraday rotation. We estimate the overall error 
in the EVPA calibration to be about 4° for all frequencies. 
We summarize the results of our EVPA calibrations in Ta- 
ble ED 

We found a systematic flux excess for all sources at 
24 GHz, presumably due to the application of a somewhat 
incorrect gain during the amplitude calibration of these 
data. Checking the total and polarized fluxes at several 
points in several sources and comparing with the results for 
the corresponding 12, 15, 22 and 43 GHz images indicated 
that the total and polarized fluxes at 24 GHz were both 
consistently about 15-20% too high. Although this will af- 
fect measurements involving the flux density at 24 GHz, it 
will not affect the inferred degrees of polarization (m) or 
X values, since these are determined using ratios. We thus 
include the 24 GHz degrees of polarization and x values in 
our analysis, but the 24 GHz measurements are not used in 
the calculations of the spectral indices. 

After the initial construction of the total intensity maps, 

1 http:/ /www. physics. purdue.edu/MOJAVE/ 

2 http:/ /www. bu.edu/blazars/VLBAproject. html 

3 http://www.astro.lsa.umich.edu/obs/radiotel/umrao.php and 
M. Aller, personal communication 
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we used the final self-calibrated visibility data to make maps 
of the distributions of the Stokes parameters Q and U at 
each frequency. Those maps were u sed to con struct the dis- 
tribution of polarized flux (p = \J Q 2 + U 2 ) and polariza- 
tion angle (x = ^arctan^), using the AIPS task COMB. 
We obtained the \ values at specified locations and their 
uncertainties using the AIPS facility IMSTAT applied to a 
3x3 pixel region (i.e., 0.15 x 0.15mas 2 , as each pixel cor- 
responds to 0.05 mas) surrounding the location of interest, 
with an estimated EVPA calibration uncertainty of 4° added 
in quadrature. 

2.2 Optical 

Nightly optical polarization observations spanned the VLBA 
observing run. Data were acquired using the SPOL spec- 
tropolarimeter (Schmidt et. al 1992a) at the 2.3m Bok tele- 
scope of Steward Observatory. On various nights, the in- 
strument was configured for spectropolarimetry using a 600 
line/mm diffraction grating. Data acquisition and reduction 
closely follow those described by Smith et. al (2003) and 
Gabuzda et. al (2006). The spectropolarimetric observations 
are averaged over the R filter bandpass for direct compar- 
ison to the imaging polarimetry. The interstellar polariza- 
tion standard stars BD+59o389, BD+64ol06, HD 236633, 
HD 245310, and Hiltner 960 (Schmidt et al. 1992b) were 
observed to calibrate Xopt- 

Optical measurements for 1633+382 are publicly avail- 
able for the Fermi progranQ. For some of the sources ob- 
served in radio, scans ranged all over the 24 hour observa- 
tion. This was the case for 0256+075, 0420-014, 0906+430 
and 1633+382. In these cases, an average of the optical 
values corresponding to days November 1 and 2, 2008 was 
taken. Both the optical and radio polarization measurements 
were corrected for statistical biases (Wardle & Kronberg 
1974). 



3 RESULTS 

3.1 Intensity and Polarization of the Cores 

The total intensity, polarization and spectral-index maps for 
all eight sources observed are presented by Algaba (2010) 
and are also available by contacting J. C. Algaba. Total and 
polarized intensity properties of the VLBI core components 
are shown in Table 2. Column [1] shows the name of the 
source, column [2] the frequency, column [3] the peak total 
intensity and column [4] the degree of polarization obtained 
by applying IMSTAT in AIPS over a 3 x 3 pixel area sur- 
rounding the intensity peak. 

We find a systematic frequency dependence in 7 out of 
8 cases, where the optical m is higher than the measured 
m at radio frequencies. However, in the radio band, there 
is no universal behaviour for the frequency dependence of 
the core m. There is some tendency for m to increase with 
increasing radio-frequency for 0133+476, 0745+241, and 
0906+430, while 0256+075 shows the opposite tendency, 
consistent with the inferred core spectral indices. There is 

4 http:/ /james. as. arizona.edu/~psmith/Fermi 



no obvious systematic radio-frequency dependence of m for 
the remaining four objects. 

All the sources have a bright core with a jet consisting 
of several components, with the innermost jet component 
not well resolved from the core at our lowest observed fre- 
quencies. The polarization angle varies with frequency in 
the core, suggesting Faraday rotation (discussed further be- 
low), while the polarization angles are normally fairly well 
aligned in the jet, indicating only modest Faraday rotation. 
This agrees with our expectation that both the magnetic 
field and the electron density decrease with distance from 
the core. 



3.2 Core Spectral Indices 

It is of interest to derive the spectral indices for the cores, 
to determine whether the core regions are predominantly 
optically thin or optically thick. Unfortunately, during the 
process of VLBI calibration and imaging, information about 
the absolute position of the source is lost. In practice this 
means that, prior to any comparison of maps obtained at 
two different frequencies, we must align the maps. This can 
be done in several ways: using astrometric phase-referenced 
observations to compare the position of the target source 
relative to a reference point-like source (see e.g., Pradel et. 
al 2006) or cross-correlating the optically thin regions of the 
jet at pairs of frequencies (see e.g., Croke & Gabuzda 2008). 
We use here the latter approach. 

For the construction of the spectral-index maps, we 
made versions of the final images using the final calibrated 
data but convolving the images at all frequencies with the 
beam obtained for the 15 GHz map. This beam was cho- 
sen as a compromise between the lower (12 GHz) and 
higher (22, 24, 43 GHz) resolution maps; although we are 
super-resolving slightly at 12 GHz, this is only by a mod- 
est amount. We then used the task COMB to obtain the 
spectral-index maps after aligning the images using the pro- 
gram developed by Croke and Gabuzda (2008). 

Table 3 contains the core shifts at each frequency rela- 
tive to 43 GHz in milliarcseconds. In each case, the derived 
core shifts are roughly opposite to the innermost jet direc- 
tions, as expected (the position of the "core" moves down 
the jet at lower frequencies) . Note that we have included the 
24 GHz data in our core-shift analysis: although the abso- 
lute scale of our 24 GHz images is inaccurate, the form of 
the image structure is correct, so that a correlation of the 
24 GHz and 43 GHz images yields a correct relative shift 
between them. The shifts between 22 and 24 GHz are small 
compared with the shifts between other pairs of frequencies, 
as would be expected, and the shifts obtained for the 24-GHz 
images are consistent with the other relative shifts obtained. 
Of course, we are not able to calculate spectral indices using 
the 24-GHz images due to the flux-scale problem described 
in §2.1. 

The core spectral indices (defined as S oc v +a ) derived 
after aligning the maps at pairs of frequencies are shown in 
Table 4. The frequencies of the two maps taken to obtain 
the data are shown in the second column. The spectral in- 
dex for the core is shown in the third column. Most of the 
source cores seem to be appreciably optically thin; i.e., their 
spectral indices are negative. The quasar 1633+382 shows 
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Table 2. Total Intensity and Degree of Polarization 



Source 


Frequency (GHz) 


T 1 / T \ 

ipcak (mjy) 


m (7o) 


I 1 ] 


[2] 


[3] 


r a i 
[4] 


0133+476 


12 


3180 ± 20 


1.1 + 0.2 




15 


3060 ± 20 


1.5 ± 0.1 




22 


2990 ± 20 


2.6 ± 0.1 




24 




2.5 ± 0.1 




43 


2370 ± 20 


2.5 ± 0.1 




Optical 




3.8 ± 0.1 


0256+075 


12 


316 + 6 


5.4 + 0.3 




15 


328 ± 4 


4.6 ± 0.2 




22 


378 ± 4 


2.6 ± 0.1 




24 




3.1 ± 0.1 




43 


389 ± 2 


1.0 ± 0.1 




Optical 




20.5 ± 0.1 


0420-014 


12 


3710 ± 20 


2.2 + 0.1 




15 


3580 ± 40 


2.4 ± 0.1 




22 


3240 ± 20 


2.4 ± 0.1 




24 




2.7 ± 0.1 




43 


2270 ± 40 


1.5 ± 0.1 




Optical 




3.9 ± 0.1 


0745+241 


12 


1390 ± 20 


1.3 + 0.1 




15 


1330 ± 10 


1.8 ± 0.1 




22 


1310 ± 9 


2.5 ± 0.1 




24 




2.6 ± 0.1 




43 


1067 ± 8 


3.5 ± 0.1 




Optical 




6.7 ± 0.1 


0906+430 


12 


694 + 6 


0.2 + 0.1 




15 


622 ± 4 


0.3 ± 0.1 




22 


590 ± 4 


1.9 ± 0.1 




24 


- 


2.1 ± 0.1 




43 


439 ± 4 


1.2 ± 0.1 




Optical 




4.5 ± 0.2 


1633+382 


12 


2110 + 20 


1.2 + 0.1 




15 


1990 ± 10 


0.4 ± 0.1 




22 


2020 ± 10 


0.2 ± 0.1 




24 




0.5 ± 0.1 




43 


1920 ± 20 


1.3 ± 0.1 




Optical 


0.33 ± 4 


3.0 ± 0.1 


1823+568 


12 


1180 + 10 


7.4 + 0.3 




15 


1131 + 7 


7.6 + 0.1 




22 


1089 ± 7 


9.2 + 0.1 




24 




8.9 + 0.1 




43 


887 + 8 


5.7 + 0.1 




Optical 




22.9 + 0.1 


1954+413 


12 


1200 ± 10 


1.4 + 0.1 




15 


1161 + 9 


1.8 + 0.1 




22 


1131 + 8 


1.1 + 0.1 




24 




1.4 + 0.1 




43 


709 + 8 


0.7 + 0.1 




Optical 




1.1 + 0.2 



Table 3. Core Shifts 


Source 


Frccjucncy Pair 


Core Shift 




USOQ y\jnZ) 


Magnitude (mas) 


Direction (°) 


I 1 ] 


[91 




[41 


0133+476 


12-43 


0.10 + 0.04 


135 + 8° 




lO io 


0.07 + 0.04 


135 + 8° 




22—43 


0.04 ± 0.04 


135 + 8° 




24—43 


n nzL -t- n c\a 


100 in o 


0256+075 


12-43 


0.16 + 0.04 


26 + 5° 




15—43 


0.11 + 0.04 


26 + 5° 




22—43 


0.06 + 0.04 


26 + 5° 




24—43 


n ofi -t- n nA 


in o 


0420-014 


12-43 


0.18 + 0.04 


+ 4° 




15—43 


0.14 + 0.04 


+ 4° 




22—43 


0.08 + 0.04 


+ 8° 




24—43 


n -t- n c\a 


n -4- a° 


0745+241 


12-43 


0.08 ± 0.04 


180 ± 8° 




15—43 


0.05 + 0.04 


116 + 10° 




22—43 


0.02 + 0.04 


90 + 8° 




24—43 


n n9 -t- n nA 
u.uz in u.u^ 


on -t- £° 


0906+430 


12-43 


0.14 + 0.04 


-45 ± 8° 




15—43 


0.11 + 0.04 


-45 ± 8° 




22—43 


0.08 + 0.04 


-45 ± 16° 




24—43 


0.05 + 0.04 


-45 ± 8° 


1633+382 


12-43 


0.10 + 0.04 


135 + 8° 




15-43 


0.07 + 0.04 


135 + 8° 




22-43 


0.04 ± 0.04 


135 + 16° 




24-43 


0.04 + 0.04 


135 + 8° 


1823+568 


12-43 


0.16 + 0.04 


26 + 5° 




15-43 


0.10 + 0.04 


26 ± 10° 




22-43 


0.05 + 0.04 


26 ± 20° 




24-43 


0.03 + 0.04 


26 ± 10° 


1954+413 


12-43 


0.13 + 0.04 


135 ± 8° 




15-43 


0.10 + 0.04 


135 + 8° 




22-43 


0.06 + 0.04 


135 + 16° 




22-43 


0.02 + 0.04 


135 + 8° 



a somewhat flatter spectrum, and the BL Lac 0256+075 an 
inverted spectrum. 



4 SEARCH FOR CORRELATIONS IN X 

We obtained the \ values of the core regions by applying 
the AIPS task IMSTAT to corresponding 3 x 3-pixel areas 
at all 5 frequencies. For consistency, we checked other nearby 
regions for each source to make sure the x values we were 
obtaining were stable around the central measured region. 
We then fitted the obtained x values with a linear A 2 law to 
obtain the Faraday rotation measure (RM) and the intrinsic 
Xo- X ~ Xo + RMA 2 . In contrast to our previous observations 
(Paper I), we now have five frequencies, two of them (22 
and 24 GHz) close enough so that we should be able to 
unambiguously determine if nn rotations are needed for the 
lower-frequency x values. 

We essentially fixed the 43 GHz value, X43, and allowed 
it to be rotated by no more than ±7t. Next, we used the 
difference between the 22 and 24 GHz x's as a guide about 
both the sign and the amount of the rotation measure. That 
is, if X22 and X24 are appreciably different, this suggests the 
presence of relatively high Faraday rotation, which must be 
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Table 4. Core Spectral Indices 



oonrcc 


x 1 rGQU.cn cy .r eiir 


opectTeil Index 




Used (GHz) 




[i] 


[2] 


[3] 




19 1 ^ 


o i s -i- o 




15-22 


-0.07 + 0.08 




22-43 


-0.34 + 0.10 


UZOO-pU i 


19 1 ^ 


in Ad A- 90 
-hu.^tu nz U.ZU 




15-22 


+0.53 + 0.14 




22-43 


+0.16 + 0.20 


0/190 01 A 
U^tZU-Ul^t 


19 1 ^ 

1Z— 10 


o 9n -i-oio 




15-22 


-0.27 + 0.11 




22-43 


-0.64 + 0.08 




19 1 ^ 
iz— 10 


O 1 Q -U O 19 

— u. iy nz u.iz 




15-22 


-0.04 + 0.10 




22-43 


-0.29 + 0.08 


uyuo-r^ou 


19 1 ^ 

IZ - 10 


O A(\ -I-OIO 




15-22 


-0.14 + 0.09 




22-43 


-0.44 + 0.09 


1633+382 


12-15 


-0.25 + 0.09 




15-22 


+0.04 + 0.08 




22-43 


-0.07 + 0.07 


1823+568 


12-15 


-0.20 + 0.11 




15-22 


-0.08 + 0.08 




22-43 


-0.30 + 0.05 


1954+413 


12-15 


-0.15 + 0.10 




15-22 


-0.09 + 0.10 




22-43 


-0.68 + 0.07 



taken into account when fitting the lower-frequency x' s - We 
then apply an appropriate number of nir rotations to Xis and 
X12 to obtain an adequate A 2 fit. In some cases, when there 
was evidence that different frequencies were probing appre- 
ciably different regions, we excluded individual end points 
from the fits. When more than one comparably good fit was 
possible, we gave preference to fits that excluded fewer data 
points and corresponded to lower RMs. As each of the cores 
is in the same optical-depth regime at all the observed fre- 
quencies, given the spectral indices (Table 4), no n/2 rota- 
tions to take into account optical-depth transitions between 
frequencies should be needed. Our results are summarized 
in Table [S] Plots of the fits given in the table are shown in 
Fig. [T] The \ values for all 5 frequencies are plotted together 
with Xopt (leftmost point in the figure) for comparison. 

Overall, it worked well to use the 22 and 24 GHz data as 
a guide to check for possible high rotation measures. The dif- 
ference between these two x values is equal to zero within la 
for 0256+075, 0420-014, 0745+241 and 1823+568; the two 
X values differ by 1.5a for 0133+476 and 2.5cr for 1945+513. 
The X22 and X24 values for 0906+430 and 1633+382 differ 
by more than 3<r. Thus, the sources that are most likely to 
display high Faraday rotation in their cores are 0906+430, 
1633+382 and 1954+513. 

In some cases, the collected \ values for the five fre- 
quencies did not lead to a single, clearly best A 2 fit. Some 
alternatives to the fits are summarized in Table [6] whose 



structure is essentially the same as Table [5] We briefly dis- 
cuss our results for each of the eight objects below. 

0133+476. The values for X43, X24 and X22 agree well 
with a linear A 2 dependence, indicating only a modest rota- 
tion measure in the core region. Further, it is possible to ob- 
tain a fit that is almost as good with a similar RM by includ- 
ing X12 together with the high-frequency points; however, 
X15 disagrees with this fit by more than 100°. We unsuc- 
cessfully tried various possibilities to see whether there was 
a high-rotation-measure fit that would provide an equally 
good fit including all 5 x measurements. We thus propose 
that the 12 and 15 GHz x values are probing a different re- 
gion in the source than the higher-frequency x values, and 
have used only the three highest radio frequencies for our 
preferred A 2 fit; the resulting Ax value is 36°, and the rms 
of the fit is less than 2° . 

0256+075. We are able to obtain an acceptable A 2 fit 
for all the points except for X43- This suggests we are probing 
a significantly different region at 43 GHz. Thus, we did not 
include 43 GHz in our final fit; the resulting Ax value is 8°, 
and the rms of the fit is about 3°. 

0420—014. There is no indication of high Faraday rota- 
tion, and we adopted the straightforward RM fit in Table [SJ 
based on all 5 frequencies; the resulting Ax value is 2°, and 
the rms of the fit is about 4°. 

0745+241. There is no indication of high Faraday rota- 
tion, and we adopted the straightforward RM fit in Table [5j 
based on all 5 x values; the resulting Ax value is 41°, and 
the rms of the fit is about 4°. 

0906+430. A good A 2 fit is not possible using all 5 x 
values, even applying reasonable rnr rotations to the lower 
frequencies. An acceptable fit is obtained if we exclude X43> 
but the rms remains fairly high, about 6°. The implied A^ 
value in this case is 2° . The three highest-frequency x values 
are very well aligned, suggesting that these higher frequen- 
cies probe a different region than the lower two frequencies, 
and so should be fit on their own. In this way, the resulting 
Ax value is 86° , and the rms of the fit is about 2° . Overall, 
it is difficult to be sure which of these should be preferred, 
and both options are given in Table [5] and shown in Fig. [T] 

1633+382. Inspection of Table[5]makes it immediately 
clear that it is not possible to obtain an acceptable A 2 fit 
with the nominal measured x values. Further, the large dif- 
ference between X22 and X24 suggests a fairly large core RM. 
Applying rotations of +ir to the Xi5 an d +3-7T to X12 yields 
a good fit to all the frequencies except for 15 GHz; at the 
same time, the parameters for the fits including and exclud- 
ing X15 are fairly similar. In addition, these fits are close to 
those for the fit that is obtained using only the x values at 
the three highest frequencies. We adopt the three-frequency 
fit as our preferred fit, assuming that the increased rms ob- 
tained when all five frequencies are used reflects the fact 
that the lower-frequency data are probing somewhat differ- 
ent regions within the source. The resulting Ax value is 33° , 
and the rms of the fit is about 3°. 

1823+568. There appears to be a clear misalignment 
between X43 and the rest of the frequencies. An excellent 
A 2 fit is obtained using the remaining four x values. This 
suggests that we are beginning to sample a significantly dif- 
ferent region at 43 GHz. Thus, we did not include 43 GHz 
measurement when fitting the x's for this core. The resulting 
Ax value is 60°, and the rms on the fit is less than 1°. 
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Figure 1. Rotation measures for radio wavelengths and optical \'s for comparison. 



1954+513. Xi'i ues substantially above the A 2 fit ob- 
tained using all five \ values. If we obtain a fit excluding this 
point, the resulting RM differs only slightly, but the RMS 
decreases by about a factor of four. We accordingly adopt 
the fit without X43- The resulting Ax value is 80°, and the 
rms of the fit is about 3°. It is intriguing that Xopt, X43 an d 
X24 are well aligned (dashed in Fig. 1), but this may be a co- 
incidence, since \22 is offset somewhat from this line. If it is 



not a coincidence, then the excellent alignment of Xopt, X43 
and X24 obviously implies Ax ~ 0, with RM = -22300+600 
rad/m 2 . 



Table 5. Measured x values and results for x vs. ^ 2 nts 



Source 


Xl2GHz(°) 


Xl5GHz(°) 


X22GHz(°) 


X24,GHz(°) 


X43GHz(°) 


X0(°) 


RM (rad/m 2 ) 


RMS 


Xo P t(°) 


Ax(°) 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


0133+476+ 


-43 ± 4 


-67 ± 4 


81 + 4 


90 + 4 


103 ± 4 


111 ± 4 


—2500 ± 500 


1.5 


75 + 1 


36 


0256+075* 


48 + 4 


19 + 4 


8 + 4 


7 + 4 


33 + 5 


—9 + 5 


1530 ± 200 


2.7 


— 17 ± 1 


8 


0420-014 


64 + 4 


42 + 4 


35 + 4 


41 + 4 


34 + 5 


29 + 5 


880 ± 220 


4.1 


27 + 1 


2 


0745+241 


133 + 4 


103 + 4 


61 + 4 


63 + 4 


56 + 4 


42 + 5 


2550 ± 270 


4.4 


83 + 1 


41 


0906+430* 


9 + 5 


65 + 7 


87 + 5 


105 + 4 


184 + 4 


129 + 9 


-3200 ± 500 


5.9 


131 ± 3 


2 


0906+430* 


9 + 5 


65 + 7 


87 + 5 


105 + 4 


184 + 4 


37 + 6 


-12100 ± 700 


2.4 


131 ± 3 


86 


1633+382* 


78 + 4 


27 + 4 


77 + 4 


42 + 6 


-100 + 21 


20 + 8 


22040 ± 960 


3.3 


-13 ± 1 


33 


1823+568* 


23 + 4 


5 + 4 


-9 + 4 


-10 + 4 


+ 4 


-22 + 1 


1250 ± 30 


0.3 


38 + 1 


60 


1954+513* 


-72 + 4 


-60 + 4 


-83 + 4 


-68 + 4 


-110 + 4 


86 ± 16 


-13800 ± 300 


3.2 


146 ± 5 


80 



Columns show [1] the source name, [2]-[6] nominal x values in degrees, [7] xo m degrees, [8] RM in rad/m 2 , [9] rms of the A 2 fit, 
[10] optical polarization angle and [11] Ax = \Xopt — Xo\- 

* XazGH z wa s not used in the fit. See text for details. 

* X\2GHz and Xi5GHz were not used in the fit. See text for details. 



Table 6. Alternative possibilities for RM fits 



Source 


Case 


Xl2GHz(°) 


Xl5GHz(°) 


X22GHz(°) 


X24GHz(°) 


X43GHz(°) 


X(.(°) 


RM (rad/m 2 ) 


RMS 


Ax(°) 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


0133+476 


A 


-43 + 4 


-67 + 4 


81 + 4 


90 + 4 


103 + 4 


120 ± 35 


-5400 ± 1800 


27.9 


45 




B 


-223 + 4 


-67 + 4 


81 + 4 


90 + 4 


103 + 4 


168 ± 22 


-10700 ± 1100 


18.7 


93 




C 


677 + 4 


473 + 4 


261 + 4 


270 + 4 


103 + 4 


84 + 22 


17100 ± 1100 


20.4 


9 




D 


-223 + 4 


-67 + 4 


81 + 4 


90 + 4 




23 + 9 


-12000 + 400 


6.0 


128 




E 


-43 + 4 




81 + 4 


90 + 4 


103 + 4 


126 + 7 


-4700 ± 400 


6.1 


51 




F* 






81 + 4 


90 + 4 


103 + 4 


111 + 4 


-2500 ± 500 


1.5 


36 


0906+430 


A 


9 + 5 


65 + 7 


87 + 5 


105 + 4 


184 + 4 


167 ± 22 


-5000 ± 1400 


19.8 


36 




13 


-351 ± 5 


-115 + 7 


87 + 5 


105 + 4 


184 + 4 


67 + 13 


-16500 + 800 


11.4 


64 




C* 


9 + 5 


65 + 7 


87 + 5 


105 + 4 




129 + 9 


-3200 ± 500 


5.9 


2 




D* 






87 + 5 


105 + 4 


184 + 4 


37 + 6 


-12100 ± 700 


2.4 


86 


1633+382 


A 


78 + 4 


27 + 4 


77 + 4 


42 + 5 


-100 + 21 


26 + 27 


22460 ± 1350 


27 


39 




B 


78 + 4 


27 + 4 


77 + 4 


42 + 5 




-48 ± 14 


28700 ± 650 


9.7 


35 




C* 






77 + 4 


42 + 5 


-100 + 21 


20 + 8 


22040 ± 960 


3.3 


33 


1954+513 


A 


-72 + 4 


-60 + 4 


-83 + 4 


-68 + 4 


-110 + 4 


86 + 16 


-14800 + 800 


13.3 


60 




B* 


-72 + 4 


-60 + 4 


-83 + 4 


-68 + 4 




66 + 6 


-13800 + 300 


3.2 


80 




C 


-792 ± 4 


-420 ± 4 


-83 + 4 


-68 + 4 


70 + 4 


174 ± 20 


-27000 ± 1000 


15.6 


25 



* Preferred fit. 

* Both fits comparable. 
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5 DISCUSSION 

5.1 The Optical and VLBI Core X Values 

The five-frequency data have partially resolved the issues 
we had with previous three-frequency data: i) we can essen- 
tially rule out internal Faraday rotation as the cause of any 
misalignments between Xopt and xo and ii) we now have a 
better basis to determine cases of high Faraday rotation in 
the VLBA cores. If any of the cores were undergoing internal 
Faraday rotation, both an oscillation of \ m the observed 
frequency range and a characteristic frequency dependence 
of the degree of polarization would be expected, but are 
not observed (see Paper I). Comparison of X22 and X24, has 
proved to be a good (although not perfect) tool to identify 
possible high Faraday rotation. 

Although we are able to obtain reasonable A 2 fits us- 
ing all 5 frequencies in some cases (0420—014, 0745+241), 
X43 deviates from the linear behaviour shown by the other 
frequencies in several other sources (0256+075, 1823+568, 
1954+513 and possibly 0906+430). This can be understood 
if the 43-GHz VLBA observations are probing a significantly 
different region than the lower-frequency observations. Sim- 
ilarly, in some sources, X43i X24 and X22 are well aligned, but 
X15 and X12 are offset from this linear relation (0133+476, 
possibly 0906+430). 

Upon inspection of Tabled it becomes clear that there 
is no obvious pattern for Ax: two sources have Ax < 10° 
(0256+075, 0420-014), one has A X > 80° (1954+513), 
and the Ax value for 0906+430 is most likely close to 90°, 
but may be close to 0°. The remaining four Ax values are 
all between 30° and 60°. There is no obvious pattern for 
Ax' = I Xopt — X43I when X43 is not aligned with other ra- 
dio frequencies. Thus, the additional information provided 
by the increased frequency coverage of our new data set has 
not led to the emergence of a clear peak in the Ax distribu- 
tion for the 8 AGN considered here. 

This has a number of implications. First, this supports 
the preliminary conclusion of Paper I that not all sources 
display aligned optical and VLBA core polarization angles. 
Second, there may be some factor differentiating the sources 
in our total sample with Ax ~ from those for which Ax 
is far from zero; e.g., they may be intrinsically different, or 
there may be some physical or observational process lead- 
ing to this difference. Third, the Ax behaviour shown by a 
particular object could in principle change with time. For 
example, D'Arcangelo et. al (2007) showed that the opti- 
cal and radio core position angles of 0420-014 followed each 
other fairly well over several days, but we have no infor- 
mation about whether this continues over long periods of 
months or years. In the following sections, we will consider 
these issues. 



5.2 Time Variation of Ax 

The observations discussed in Paper I and those presented 
here yield a subset of sources for which we have Ax mea- 
surements for two epochs; these results are summarized in 
Table 7. 

The Ax values inferred for the two epochs are roughly 
the same for 0133+576, 0745+241, 1633+382, 1823+568 and 
possibly 0906+430, whereas the Ax values for the other 



sources differ by more than about 20°. Thus, roughly half 
of these sources display a substantial change in their Ax 
estimates over time scales of several years. Comparing the 
core RM fits for the two epochs (columns [5] and [6]), we do 
not see any clear trend for systematically higher or lower RM 
values in our five- frequency data (epoch D), compared to the 
three-frequency data (epochs GA, GB, A, B, C). Overall, 
the RM values are similar to order of magnitude, with the 
exception of 1633+382, which differ by a factor of 50. This 
suggests that there were no systematic problems with our 
previous three-frequency RM fits, and that those fits were 
generally reliable despite the limited frequency coverage. 

The signs of the core rotation measures for 0133+476, 
0745+241 and 1633+382 apparently changed between 
epochs. However, this is not unprecedented. Zavala & Taylor 
(2001) found similar changes in the core rotation measures 
in 3C273 and 3C279. 

If Ax changes over time, we can consider several pos- 
sibilities. First, it might be that Xopt and xo each change 
independently, such that Ax varies randomly. It could also 
be that Ax follows some predictable path, such as an oscil- 
lation around a certain value or a rotation over time. Mon- 
itoring of Xopt and xvlba over a longer time period would 
be required to test this. A third option would be that Ax is 
correlated with some other physical properties of the source 
- for example, that Ax is close to zero only when the core is 
active, or quiescent. We consider this last possibility in the 
following sections. 



5.3 Other Physical Properties and Ax 

In principle, a characteristic difference between the optical 
and the Faraday-corrected VLBA core polarization angle 
could arise as the result of some physical properties of the 
core region. We considered several possibilities in Paper I. 
For example, we found no evidence for a correlation between 
Ax and |xo — where Qj is the inner jet direction. No 
evidence was found for a correlation between the degree of 
polarization in the core and Ax, although there were slight 
indications that Ax varies over a smaller range of values as 
the degree of optical polarization m op t increased. Finally, we 
found no link between depolarization and Ax- 

We consider two further possible relations between Ax 
and other physical characteristics of the source, namely, the 
intrinsic core RM and the core-region magnetic field. For 
completeness, we will do this for all available epochs, in- 
cluding those analyzed in Paper I. 



5.3.1 Rotation Measure 

We first look for a correlation between the rotation mea- 
sures found for our cores and Ax- For sources with known 
redshift, intrinsic rest frame core rotation measures can be 
calculated using RM int = RM b s (i + z) 2 . The results are 
given in Fig. [51 where black squares indicate BL Lac ob- 
jects; dark gray diamonds high polarized quasars (HPQs) 
and light grey stars low polarized quasars (LPQs) (for a 
complete list of the RMs, see Paper I). It is clear that no 
correlation appears. The magnitude of the rotation measure 
has no obvious relationship to the observed A\- 
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Table 7. Time evolution of core RM and Ax 



Source 


Epoch 1 


Epoch 2 


At (yr) 


RM (Ep. 1) (rad/m 2 ) 


RM (Ep. 2) (rad/m 2 ) 


A x (Ep. 1) (°) 


A x (Ep. 2) (°) 


[11 


f2l 


[31 


[41 


[51 


[61 


[71 

L'J 


[81 

L U J 


0133+476 


c 


D 


3.10 


200 ± 30 


-2500 ± 500 


59 


36 


0256+075 


A 


D 


4.01 


200 ± 820 


1530 ± 200 


60 


8 


0420-014 


B 


D 


3.64 


270 ± 200 


880 ± 220 


43 


2 


0745+241 


B 


D 


3.64 


-1780 + 320 


2550 ± 270 


51 


41 


0906+430t 


A 


D 


4.01 


-24440 ± 470 


-3200 ± 500 


8 


2 


0906+430* 


A 


D 


4.01 


-24440 ± 470 


-12100 ± 700 


8 


86 


1633+382 


A 


D 


4.01 


-570 ± 430 


22040 ± 960 


29 


10 


1823+568 


GA 


D 


6.24 


400 ± 100 


1250 ± 30 


54 


60 


1954+513 


B 


D 


3.64 


-24400 ± 640 


-13800 ± 300 


4 


60 


OJ287 


GB 


C 


2.56 




-1920 + 120 


3 


68 


3C279 


GB 


B 


2.03 


-2360 ± 120 


-1020 + 110 


3 


18 



Columns are as follows: [1] source name, [2]-[3] epochs observed (GA: August 7, 2004; GB: March 5, 2003; A: November 1, 2004; 
B: March 15, 2005; C: September 26, 2005; and D: November 2, 2008), [4] time (in years) between those two epochs, [5]-[6] derived 
RMs and [7]— [8] the two Ax values at the corresponding epochs. 

* Fit obtained using the x values at 12, 15, 22 and 24 GHz. 

* Fit obtained using the x values at 22, 24 and 43 GHz. 
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Figure 2. Intrinsic rotation measure RM compared against Ax- 
Black squares, dark gray diamonds and light gray stars indicate 
BL Lac objects, HPQs and LPQs respectively. Note that, for clar- 
ity, 0906+430, 2230+014 for epoch A and 1954+413 for epoch D 
have not been plotted as their intrinsic RM is higher than 60 000 
rad/m 2 . 



5.3.2 Core-region Magnetic Field 

Typically, due to the finite resolution of the observations, the 
VLBI "core" is in fact a combination of the true optically 
thick core (i.e., the surface with optical depth r = 1) plus 
the unresolved innermost region of the optically thin jet. Due 
to the frequency dependence of the r = 1 surface (Konigl 
1981), the position of the radio core varies as r oc v~ 1 / kr , 
with r the distance from the central engine and k r — ((3 — 
2a)m + 2 — 2)/(5 — 2a), where m and n describe the decay of 
the magnetic field (B — Bir m ) and particle number density 
(N = Nir n ) with distance from the jet base (i.e., the apex 
of the cone formed by the jet). In the case of equipartition, 
k r — 1, independent of the spectral index a, with the choices 
of m = 1, n = 2 being physically reasonable. 

We can derive the magnetic field strength from this core 
shift combining the results of Lobanov (1998) and Hirotani 
(2005) (also see O'Sullivan & Gabuzda 2009). The magnetic 
field in Gauss at 1 pc from the jet base is given by 



B 1 ~ 0.025 



1/4 



(1) 



where z is the redshift, 5 the Doppler factor, <j) the half 
opening angle of the jet, 9 the viewing angle and Q rv the 
core position offset, defined by Lobanov (1998) as 



£l rv = 4.85 x 10" 



1 + z 2 



l/k r l/k r 



l/k r 



l/k T 



(2) 



with Dl being the luminosity distance and Ar mas the core 
shift in mas found between the frequencies v\ and v^. 

We have estimated k r for the sample using an expo- 
nential of the form Ar = A x (v~ 1/kr - (43)~ 1/fcr ), where 
Ar are the core shifts given in Table 3 and 43 GHz is used 
as the reference frequency. Fits are shown in Figure and 
the corresponding k r values are given in Table [8] For all 
but two sources (0906+430 and 1823+568), k r is equal to 
1 within the la uncertainties, and for only one (1823+568) 
is k r 7^ 1 within 2a. Thus, the values obtained for k r are 
mostly compatible with equipartition in the core region. 

Once the core shifts and k r have been obtained, the core 
position offset Q, „ can be calculated. As the majority of the 
cores are found to be close to equipartition, the magnetic 
field B at 1 pc from the jet base can be derived, provided 
that there are adequate values for the Doppler factor, the 
jet opening angle, the viewing angle and the redshift. These 
parameters and the inferred results are also given in Table [8] 
In all cases, a jet opening angle c/> = 0.5° is assumed. 

Values of Bi were obtained for four sources (Table |SJ). 
These values are typically ~ 0.7 Gauss, similar to previous 
results for other radio-loud AGN (Lobanov 1998, O'Sullivan 
& Gabuzda 2009) . There is no tendency for cores displaying 
k r values far from unity (indicating that they are relatively 
far from equipartition) to stand out in terms of their Ax 
values (Figure 4). 

A check was made for any dependence of Ax on Q. rv , 
as derived from the core shifts, and on the strength of the 
magnetic fields deduced. However, there are five-frequency 
data for only eight sources, and we have sufficient infor- 
mation to derive magnetic fields for only half of these; too 
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Figure 3. Frequency dependent shift of the core position for the sources studied. Fits correspond to the spectral shift power index k r 
when 43 GHz is used as the reference frequency. 



few sources to be of any statistical value. This issue can be 
partially addressed by deriving S} rv (and hence, magnetic 
fields) for the sources studied in Paper I. However, those 
data were obtained at only three frequencies (15, 22 and 
43 GHz), which proved to be insufficient for reliable fits for 
k r - Therefore, based on our five- frequency results and the re- 
sults of O'Sullivan & Gabuzda (2009), we assume equiparti- 



tion (k r ~ 1) for the observations in Paper I; this assumption 
is likely to be incorrect only for a small minority of those 
sources. The equivalent of Table [8] for the data of Paper I 
is summarized in Table [5] where sources have been marked 
as belonging to epochs A (November 1, 2004), B (March 15, 
2005) or C (September 26, 2005). 

The magnetic fields obtained range from ~ 0.2 Gauss 
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Table 8. k r , Q. TV , and B\ 



Source 




z 


D L (Mpc) 


*(°) 


rt 


6(°)t 


5t 


Slrv (PC/Hz) 


Bi (G) 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


0133+476 


0.9 + 0.1 


0.86 


5.49 x 10 3 


0.5 


14.4 


2.5 


20.7 


(1.47 + 0.04) x lO 10 


0.3 + 0.3 


0256+075 


1.0 + 0.1 


0.89 


5.60 x 10 3 


0.5 








(2.20 + 0.20) x 10 10 




0420-014 


1.2 ± 0.4 


0.92 


5.93 x 10 3 


0.5 


11.4 


1.9 


19.9 


(2.11 ± 0.06) x 10 10 


0.8 + 0.3 


0745+241 


0.9 ± 0.1 


0.41 


2.22 x 10 3 


0.5 








(0.76 ± 0.03) x 10 10 




0906+430 


2.2 ± 1.0 


0.67 


6.83 x 10 3 


0.5 








(2.60 + 0.60) x 10 10 




1633+382 


0.9 ± 0.1 


1.81 


1.39 x 10 4 


0.5 


30.5 


2.5 


21.5 


(1.63 ± 0.04) x 10 10 


0.7 + 0.3 


1823+568 


0.6 ± 0.1 


0.66 


3.99 x 10 3 


0.5 


37.8 


5.0 


6.4 


(1.70 + 0.60) x 10 10 


0.7 + 0.2 


1954+413 


0.9 ± 0.4 


1.22 


8.53 x 10 3 


0.5 






7.4 


(1.40 + 0.30) x 10 10 





Columns are as follows: source name, spectral shift power index k r \ redshift z; luminosity distance Dl in Mpc; 
Lorentz factor T, viewing angle 8, Doppler factor <5, core shift parameter S7 r „ and the inferred magnetic field at 
1 pc from the jet base B\. 
t from Hovatta et al (2009) 



Table 9. Magnetic-Field Parameters for Sources in Paper I (assuming k r = 1) 



Source 


Epoch z 


D L (Mpc) 


rt 


0(°)t 


5t 


£trv (pc/Hz) 


Bi (G) 


[1] 


[2] [3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 



0133+476 


C 


0.86 


5337 


14.4 


2.5 


20.7 


(0.9 + 1.3) 


x 


10 10 


0.3 + 0.3 


0138-097 


c 


0.73 


4381 








(2.3 ± 1.0) 


x 


10 10 




0256+075 


A 


0.89 


5602 








(5.4 + 1.2) 


X 


10 10 




0420-014 


B 


0.92 


5770 


11.4 


1.9 


19.9 


(6.2 ± 1.1) 


X 


10 10 


1.7 + 0.3 


0745+241 


B 


0.41 


2164 








(6.2 ± 1.2) 


X 


10 10 




0804+499 


A 


1.43 


10427 


17.8 


0.2 


35.5 


(1.0 + 1.2) 


X 


10 10 


1.2 + 0.9 


0814+425 


C 


0.25 


2952 


2.7 


8.4 


4.6 


(1.0 + 1.1) 


X 


10 10 


0.4 + 0.2 


0859+470 


B 


1.46 


10455 








(1.0 + 1.1) 


X 


10 10 




0906+430 


A 


0.67 


3934 








(3.7 + 1.3) 


X 


10 10 




0953+254 


B 


0.71 


4208 






4.3 


(3.6 ± 1.0) 


X 


10 10 




1055+018 


B 


0.89 


5593 


11.1 


4.7 


12.2 










1156+295 


A 


0.73 


4332 
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t from Hovatta et al (2009) 



(OJ287, epoch C) to ~ 2.3 Gauss (2134+004, epoch A), with 
the typical magnetic field strengths being of the order of half 
a Gauss. Comparing the inferred magnetic fields for the two 
sources available for both epochs reveals that B\ is the same 
within the uncertainties for 0133+476, but differs by about a 
factor of two for 0420-014 (this difference is ~ 2.1<r) between 
our two epochs. 

Comparisons of the magnetic field strength with Ax are 
shown in Fig. [S] The correlation coefficients are found to be 
-0.17 for Q. rv and -0.38 for B\. However, the maximum ob- 
served Ax tends to decrease as the magnetic field increases. 



If this is true, this suggests that large misalignments in Ax 
could be associated in part with relatively low core magnetic 
fields. 

A plausible explanation is that as the magnetic field in- 
creases, any turbulence, random component or disruption in 
the AGN would be of less relative importance. Under those 
circumstances, it is reasonable to suppose that optical and 
radio properties might become more and more correlated 
and thus, Ax will tend to be small. Another way to look at 
this hypothesis is that the jet magnetic fields are not sum- 
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Figure 4. Spectral shift power index k r parameter for the Ax 
values obtained in the sources studied. No tendency in Ax is seen 
for these cores with k r differing from unity. 
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Figure 5. Core shift parameter Q rv (upper) and magnetic field 
at 1 pc (lower) plotted against Ax- Symbols have the same sig- 
nificance as in Figure 2. 



ciently well ordered over sufficiently large regions to ensure 
that Ax is always small. 



6 CONCLUSIONS 

Our results in Paper I indicated that the difference be- 
tween the VLBA radio core and optical polarizations of 



AGNs showed a more complex behaviour than was expected 
based on the results of Gabuzda et. al (2006). We have 
considered here simultaneous observations at 5 frequencies 
(12+15+22+24+43 GHz) and in the optical, in order to de- 
termine if effects such as internal or high external Faraday 
rotation could be affecting our radio data. 

With these new data, we have been able to resolve some 
ambiguities in the rotation measure, and confirm the pres- 
ence of relatively high Faraday rotation in some of the ob- 
served cores. However, no extremely high rotation measures 
were found, and the values derived from our five-frequency 
data are generally consistent with previous estimates based 
on three frequencies. 

The behaviour of the polarization angles can be com- 
plex. Some deviations from A 2 behaviour are found at high 
frequencies (usually for 43 GHz), restricting the frequency 
range that can effectively be used to derive the core RMs. 
It is possible that these deviations arise because different 
regions in the source or different source structures are being 
probed at different frequencies. 

Using five-frequency data to estimate the zero- 
wavelength radio core polarization angles for the eight 
sources considered here and comparing these with the si- 
multaneously measured optical polarization angles yields no 
clear peaks in the distribution of Ax; in particular, no clear 
peak at Ax ~ has emerged for the 8 AGN considered here. 
We have attempted to identify some characteristic property 
of the cores for which Ax ~ 0, but without success. We have 
extended the analysis in Paper I by examining the data cor- 
relations between Ax and the core rotation measures and 
the magnetic fields at 1 pc from the jet base. We have found 
no evidence for correlations between any of these properties 
and the observed Ax values. Thus, although the overall dis- 
tribution of Ax for all 40 sources in our sample does show 
a significant peak at Ax ~ (Paper I), it remains unclear 
what distinguishes these AGN cores from those having mis- 
aligned polarization position angles. 

Thus, Ax ~ for some AGN, consistent with co-spatial 
emission in the optical and radio, Misaligned polarization 
position angles possibly indicate that the optical and radio 
emission in those objects arises in different parts of the jet 
with different magnetic field orientations. There are cases 
where Ax remains constant for a particular source, but there 
are also counterexamples of the optical and VLBI core po- 
larization angles changing with time. We do not currently 
have sufficient information to identify the origin of this be- 
haviour, and further coordinated optical and VLBA obser- 
vations would be of interest in this connection. 

We have estimated the core-shift exponents k r for the 
eight AGN considered here. The fitted k r values for six of the 
eight AGN are equal to unity within the uncertainties, sup- 
porting the previous conclusion of O 'Sullivan & Gabuzda 
(2009) that the core regions of radio-loud AGN are often 
near equipartition. The derived core magnetic fields at 1 pc 
from the jet base for four AGNs are of the order of half 
a Gauss, similar to the values found for other sources by 
O'Sullivan & Gabuzda (2009), based on VLBA data ob- 
tained simultaneously at seven frequencies. 

We have found no evidence for clear correlations be- 
tween any of these properties and the observed Ax values, 
although there is a hint that the maximum observed Ax 
may tend to decrease as the core-region magnetic field in- 
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creases, suggesting that large misalignments in Ax could be 
more common in cores with relatively low magnetic fields. 
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